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ABSTRACT 


We study cosmological consequences of spontaneous breaking of an approximate discrete 
symmetry. The breaking leads to formation of proto-domains of false and true vacuum 
separated by domain walls of thickness determined by the mass scale of the model. The 
cosmological evolution of the walls is extremely sensitive to the magnitude of the biasing; 
several scenarios are possible, depending on the interplay between the surface tension on the 
walls and the volume pressure from the biasing. Walls may disappear almost immediately 

after they form, or may live long enough to dominate the energy density of the Universe 
• INFN, Italy. On leave of absence from Dip. di Fisica, Univ. di Roma II (“Tor Vergata”). 
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and cause power-law inflation. We obtain limits on the biasing that characterize each » 
possible scenario. 
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1. Introduction 
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It has been known for some time that spontaneous breaking of a discrete symmetry leads to 
formation of surface-like topological defects called domain walls. 1 These two-dimensional 
structures are stable if the vacuum manifold determined by the interaction potential is 
topologically non-trivial in the sense that it is defined by disconnected points related to 
discrete, degenerate ground states of the theory. Once the field settles into one of the 
possible ground states of the theory (i.e., a point in the vacuum manifold) it cannot 
continuously be transformed into another point in the vacuum manifold. The transition 
region between degenerate ground states is the domain wall. In the rest of this section we 
very briefly review these ideas. 

The appearance of walls in a theory with discrete symmetry breaking can be very easily 
understood with a simple model: Consider a real scalar field a with Lagrangian density 

£ = \i d n°) 2 ~ - °l) 2 - (i) 

The Z 2 symmetry of the Lagrangian is broken when a obtains a vacuum expectation value 
a — +cr 0 or a = —oq. The transition region between the two possible vacuum values for a 
defines the domain wall. 

Domain walls can be described by the solutions to the equation of motion obtained 
from the Lagrangian of Eq. (1), 

+ A<t(<t 2 - oj) = 0. (2) 

Imagine an infinite wall in the y — z plane at x = 0. The solution with boundary conditions 
a = —cq at x = — 00 , and a = +cr 0 at x = + 00 , is simply 

o(x) = Oo tan h(x/A) (3) 

with A = (A/2) 1 / 2 a^ _1 defined as the “thickness” of the wall. Balance between the potential 
energy that tends to make the wall thinner, and the gradient term that tends to make the 


3 


ORIGINAL PAGE IS 
OF POOR QUALITY 

wall thicker, gives rise to the finite thickness of the wall separating regions of different 
vacua. 

The stress tensor for the wall, T M „ = —d^ad^a — Cvnw> is 

T^v = cosh~ 4 (x/A)diag(l,0, 1, 1). (4) 

From we can define the surface tension of the wall, i.e. the energy per unit area in 
the rest frame of the wall, 

H = J T° 0 dx = J [^(V<r) 2 + ^(<r 2 - ol) 2 ]dx, (5) 

which is given by 

„ . H^a>/V 0 3 . (6) 

Since T. D. Lee’s model for CP-violation of the early seventies, 2 there have been many 
models where the presence of a real scalar field with quartic self-interaction leads to the 
spontaneous breaking of a discrete symmetry. For example, axion models which have been 
proposed as a solution to the strong CP problem in QCD have a spontaneously broken Zjq 
that generates discrete degenerate vacua. 3 We mention these two models because they were 
particularly important to the development of research on the role of topological defects in 
cosmology. 

The first effort to understand the cosmological effects of spontaneous breaking of a 
discrete symmetry was the work by Zel’dovich, Kobzarev, and Okun. 4 They outlined most 
of the details for the cosmological evolution of a Universe filled with walls in the case of 
degenerate vacua. Their conclusions were quite dramatic: The wall contribution to the 
energy density of the Universe will quickly overpower the radiation contribution, causing 
a period of power-law inflation (of course, their work pre-dated the proposal of infla- 
tion) with the scale factor going as R(t) ~ t 2 . An expansion rate that fast would leave 
less time for galaxy formation and change the production rates during nucleosynthesis. 
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Also, the presence of a wall in the observable Universe would cause distortions in the 
cosmic microwave background that would violate the present limits on its homogeneity 
and isotropy, unless the wall’s gravitational redshift is safely below the detectable limit of 
6T/T < 0(10 -4 ). These limits imply an upper bound on the surface tension fx < 0(10 -1 ) 
g cm -2 or A 1 / 6 oq < 0(10 -2 ) Gev. 4 Thus, unless these limits could be satisfied (or some- 
how walls disappear early enough in the cosmological expansion), models with discrete 
symmetry breaking are ruled out by cosmology. 

As the above limit on the scale for symmetry breaking seems very restrictive from the 
particle physics side, attention turned to the other way to accommodate walls, namely to 
have them exist for only a brief period. In fact, Zel’dovich, Kobzarev and Okun suggested 
that if the discrete symmetry was not exact, the energy difference between the two vacua 
would cause the false vacuum regions to disappear, possibly before the walls could dominate 
the energy density of the Universe. This same idea was mentioned in the works of Kibble, 5 
Vilenkin, 6 and Sikivie. 7 A small bias favouring one vacuum over the other (or others) 
can indeed make the walls go away. In Refs. (6-7) a lower bound for the asymmetry 
was obtained based upon the requirement that the walls should disappear before they 
dominate the cosmological evolution. Nevertheless, as we will point out in this work, the 
dynamics of the walls subject to both the usual tension force due to curvature and to the 
volume pressure due to the biasing is much richer than it has formerly been appreciated. 
The two forces will clearly compete with each other, with the resulting wall dynamics 
extremely sensitive to the magnitude of the bias. We will study in detail the cosmological 
consequences of the breaking of an approximate discrete symmetry (a “quasi-symmetry”). 

Different scenarios we consider may also be of relevance in the case of hybrid topological 
structures such as walls bounded by strings, 1 or in the recently developed scenario for 
formation of non-topological solitons (hereafter NTSs). 8 The case of NTSs is particularly 
interesting; the real scalar field that gives rise to the walls couples to a complex scalar 
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field (or a fermion field) that carries a conserved global charge in such a way as possibly to 
allow for a small violation of the discrete symmetry (in fact it is demanded in the model 
of Ref. 8). In the false vacuum (the NTS interior) the complex scalar field is massless 
while in the true vacuum it is massive. The two regions are separated by a domain wall of 
thickness set by the mass of the real scalar field. 

Due to the many possible situations of interest, we will not analyse any particular 
model. We will try, instead, to keep our results as general as possible. In Fig. 1 we show 
the general shape of the potential that we will consider. More details about the general 
model are given in Section 3. 

The paper is organized as follows. In Section 2 we develop the kinematics of walls by 
considering a perfect gas of walls moving with an average velocity u; in a box of volume V. 
We will obtain a velocity-dependent equation of state for the wall gas. Next, by assuming 
that the energy density of the Universe is dominated by the wall gas, we will obtain the 
velocity-dependent law of cosmological expansion. In particular, we will show that the 
maximum wall velocity that results in power-law inflation is ti> = 1/ y/3, the sound speed 
in a relativistic gas. This section can be read quite independently from the rest of the 
paper. In Section 3 we will describe the formation of walls in a primordial phase transition 
in the presence of biasing. The dynamics of walls will be studied in Sections 4 and 5, 
assuming, in section 4, that the volume pressure acts before the tension force can stretch 
the walls to the horizon scale, and, in Section 5, assuming the biasing is small enough to 
allow the tension force to straighten walls to the horizon scale. In the second case, walls 
may dominate the energy density of the Universe, and cause a power-law inflation. Finally, 
in Section 6 we will review our results and discuss possible future directions for further 
work. 
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2. Kinematics of Walls and Evolution of a Wall-Dominated Universe 


In this section we obtain the equation of state for a perfect gas of walls moving with average 
velocity tv inside a box of volume V > , where A was defined in the introduction as the 

wall’s thickness. The perfect gas assumption means that we will neglect possible dissipative 
effects that may come from interactions between walls. 

Consider N walls in the y — z plane, moving in the x direction with average velocity 
w. Each wall’s position is described by the point x % , with the index i running from 1 to 
N. If there are many walls inside the box such that the average wall separation is much 
smaller than the linear dimensions of the box, we can write an average energy-momentum 
tensor for the N walls as 

/ dxf(x)T t ti/(x) , . 

" Jdxf(x) ’ [l) 

where /(x) is the distribution function for the walls, / f(x)dx = N. If (L) is the average 

wall separation, we can approximate {T^ v ) as 

(T„„) = (jy / T^dx = (8) 

The tensor can be understood as the average energy-momentum surface density of the 

walls. For example, Sqq = /z, where \i is the surface tension of the walls defined in Eq. (5). 
It should not be confused with the spacetime index //. Thus, the wall energy-density is 
given by p w = { T 00 ) = n/{L). 



The walls will be moving with average velocity tv in the +x direction with respect to 
an observer at rest with the box. Accordingly, the tangential components of v (i.e. S 22 
and S33) will not be affected by the walls’ motion. 4 The same cannot be said of the other 


components. Upon Lorentz transformation, Sfu, for an observer at rest is 


SfiV — 


7 2 H 
y 2 fitv 
0 
0 


7 2 Htv 

7 2 /zu; 2 

0 

0 


0 

0 

0 





(9) 
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where 7 = (1 — to 2 ) -1 / 2 as usual. Of course, there will be walls moving in the —x direction 
as well. Once we average the two directions, the off-diagonal terms disappear. In the most 
general case, the same procedure must be repeated for the ±y and for the ±z directions. 
The final result for the average energy-momentum tensor for a gas of walls moving with 
velocity w with respect to an observer at rest is then 


^ - sfz> 


0 

0 


\ 


0 

0 


( 10 ) 


/ 3-7 2 0 

0 ( w 2 ') 2 - 2) 

0 

V 0 

Defining the energy density as pw = (Too) and the pressure as pw = (Tj t ), the equation 
of state for the wall gas becomes 


0 
0 

(tu 2 7 2 - 2) 0 

0 (tu 2 ^ 2 - 2) J 


(w 2 7 2 — 2) 

pw = — ^2 — p W' 


( 11 ) 


We point out three cases of interest: For w relativistic (tu ~ 1, 7 » 1), pw = 
PW 1 3, and the walls behave like a relativistic gas; for w non-relativistic (w ~ 0, 7 ~ 1), 
PW = —2pw/3, which is the well-known result for static walls; and finally, for w zt 0.82 
(tu 2 7 2 = 2), pw = 0, and the walls behave as pressureless matter. Thus, for w > 0.82 the 
gas of walls will contribute positive pressure to the energy-momentum tensor. 


Once we have the equation of state for the wall gas, it is natural to ask how the evolution 
of a Universe filled with walls scales with the velocity of the walls. For simplicity we take the 
flat Robertson- Walker metric as the spacetime metric, ds 2 = —dt 2 +R(t) 2 (dx 2 +dy 2 +dz 2 ). 
With this metric, and the equation of state given by Eq. (11), the energy-momentum 
conservation equation, d(pwR z ) = -pwd{R 3 ), can be integrated to give 

(12) 

The equation that governs cosmological evolution, ( R/R ) 2 = (8nG/3)pw , can be inte- 
grated using Eq. (12) to give 

R{t) ~ t 2 / 3 (“ +1 ). 
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The three cases of interest mentioned above will clearly correspond to the well-known 
solutions R ~ i 1 / 2 for a = 1/3; R ~ t 2 for a = —2/3; and R ~ $ 2 / 3 for a = 0. 

Notice that for walls moving slowly enough, there may be a power-law inflation when- 
ever the scale factor expands faster than t. It is easy to verify from Eq. (13) that the upper 
limit for the wall speed in order to have inflation is tv = l/y/S, the speed of sound for a 
relativistic gas ! 

Before we move on to discuss the formation of walls, we would like to comment on 
the motion of walls in the presence of two fluids that may have (and do in most cases of 
interest) different “chemistry”. In the case only the vacuum contributes to the energy- 
density of the Universe, the volume pressure accelerates the walls to c (c is the velocity 
of light ‘in vacuo’). The important point is that walls that are accelerated due to the 
difference in vacuum energy between the two sides may achieve a terminal velocity smaller 
than c in case there is radiation or non-relativistic particles on both sides. As discussed in 
the work of Steinhart, 9 this velocity is, however, larger than l/y/3 in most cases. 

Steinhart 9 has considered the motion of accelerated plane walls in the presence of two 
fluids, relativistic or not, at different temperatures on both sides of the wall. Only when 
the temperature of the fluid in the false vacuum is zero will the walls be accelerated to 
c. Using the language of “detonation waves” , the wall can be understood as a wave front 
separating the burnt fluid, i.e. the true vacuum at the back of the wall, and the unburnt 
fluid, i.e. the false vacuum in front of the wall. The equivalent of the chemical energy 
stored in the unburnt fluid in a usual detonation is then the vacuum energy stored in the 
false vacuum. As the wave front propagates, it converts false vacuum into true vacuum, 
with the energy difference being used to further accelerate the wall and to heat up the 
fluid left behind. The wall drags with it the fluid just burnt, which will not be at rest 
with respect to the Universe, contrary to the thermalised fluid further away from the wall. 
Because the velocity of the wall with respect to the fluid immediately in front of it turns 
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out to be larger than the speed of sound in the unburnt fluid, 9 a point in the unburnt fluid 
cannot be perturbed by the wall until the wall passes by it. Thus, the unburnt fluid is at 
rest with respect to the cosmic rest frame. 

Steinhart gives the velocity of the wall with respect to the unburnt fluid in contact 
with it, which is precisely the wall velocity u that appeared in the equations above. The 
terminal velocity of the wall is in general larger or equal to l/>/3; this is the smaller possible 
velocity in the case that both fluids are relativistic. When the burnt fluid is relativistic 
and the unburnt fluid is non-relativistic the minimum terminal velocity is even larger, 
w = y/Z/2. In the opposite case of a non-relativistic burnt fluid and a relativistic unburnt 
fluid the walls may move in the opposite direction, towards the true minimum side, if the 
vacuum energy cannot compensate the thermal pressure of the relativistic fluid in the false 
vacuum. This possibility has been raised in Ref. 8, where it was shown that for reasonable 
values of the couplings the thermal pressure is dominated by the vacuum pressure, causing 
the wall to move towards the false vacuum. Recall also that as the Universe expands 
the thermal pressure becomes progressively weaker, with the walls eventually reaching 
relativistic speeds. 
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We now study the spontaneous breaking of an approximate discrete symmetry in the early 
Universe. First, some comments about the potential shown in Fig. 1. For the purpose 
of this paper, the exact way in which the asymmetry appears is immaterial, so long as 
the energy difference between the two vacua can be written as A = ecq, where e can 
be a function of coupling constants and temperature. One example was given in Ref. 8. 
In the limit e — ► 0 we recover the exact degeneracy of Eq. (1). The parameter e is the 
dimensionless asymmetry parameter. 

Of course, the potential of Fig. 1 is the zero-temperature potential. It is only valid 
for temperatures safely below the critical temperature for the phase transition, T c , which 
can be calculated for a particular model. For the hot big-bang model of the Universe 
it is reasonable to assume that temperatures higher than T c were achieved sufficiently 
early in the cosmological expansion. In this case one must use the temperature-corrected 
potential. Let us consider the e = 0 symmetric potential written in Eq. (1). For T T c 
it is well known that the discrete symmetry will be restored and the potential will have 
a minimum at a = 0. The details of the phase transition for the symmetric case and the 
consequent formation of walls have been worked out by Kibble. 5 Here we repeat some of 
his arguments so that we can compare them to the asymmetric case. Most of the results 
for am approximate symmetry have been derived in Ref. 8. 

First, it is not difficult to show that for the potential in Eq. (1) the critical temperature 
is T e = 2oo- At zero temperature the two possible ground states for a axe given by a = 
±<7q. As the temperature approaches T c , thermal fluctuations in the a field become large, 
with regions rapidly (compared to the cosmological time scale) interconverting between 
the two possible zero temperature values for a. At high temperatures there is enough 
thermal energy in the system for the fluctuations to “jump over” the potential barrier. 
(At T = 0, the height of the potential barrier is given by Um — (A/4)<7q.) The typical 
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volume of a fluctuation region is given by = £ 3 , where £ is the correlation length, 
given approximately by the inverse temperature-dependent mass of the a field, £(T) _1 = 
A 1 / 2 <tq (1 — T 2 /T 2 ) 1 / 2 . At T = 0, £ = A -1 / 2 Oq 1 . Below T c the transition rate between the 
two vacua is proportional to exp {-Fm/T), where Fm is the free energy of the fluctuation, 
Fm = Um xTf. As the temperature decreases, the barrier between the two vacua increases 
while the thermal energy that drives the fluctuations in the a field decreases. One can 
define the Ginzburg temperature, Tq, as being the temperature below which fluctuations 
over the barrier will be exponentially suppressed and the population in the two vacua will 
be frozen. If one ignores the expansion of the Universe (in a more complete treatment one 
should obtain Tq by comparing the thermal fluctuation rate, Ty, to the expansion rate 
of the Universe, H : Ty ~ H, at T = Tq), the Ginzburg temperature can be estimated 
to be Tq = V^Um- 5 In the case under consideration, Tq = ao(A + 1/4) -1 / 2 . Due to the 
perfect symmetry between the two vacua, the probability for a fluctuation to end up at 
+< tq or —oq clearly is 50%. Space will then be divided into cells of volume approximately 
given by £ 3 , with walls separating cells of positive and negative vacua. Of course, between 
cells of the same vacuum there will not be a domain wall. What will then be the general 
structure of space as the phase transition is completed? This question has been answered 
by studies of percolation in large lattices. 10 It has been shown that if the probability for, 
say a plus-cell, is bigger than a certain value p c , an infinite (in an infinite lattice) plus- 
cluster appears, while if the probability is smaller than p c only finite plus-clusters appear. 
The value for p c , the percolation threshold, varies with the type of lattice considered, but 
in all 3-dimensional lattices is smaller than 50%. In the perfectly symmetric case both 
vacua have probabilities above p c ; space will be permeated by an infinite wall of very 
complicated topology dividing regions of plus and minus vacuum. There will also be small 
clusters but these will be exponentially suppressed. 

The introduction of a small bias will change slightly the arguments above. For T » T c 
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the potential will still have a parabolic shape but its minimum will not be exactly at a = 0. 
Its location will depend on the details of the model under consideration. Also, quantities 
such as fi, T c , Tq and £ will have corrections proportional to the asymmetry parameter 
e. However, we will assume that such corrections will not be very important due to the 
smallness of e, and that it is a legitimate approximation to take for these quantities the 
same values as in the e = 0 case. The same cannot be said about the relative probability 
of having a fluctuation end up in a plus cell (denoted p+) or in a minus cell (denoted 
p-). The values for p+ and p_ are very sensitive to the energy difference between the 
two vacua, A = euq. As A increases, the false vacuum, with larger free energy, becomes 
progressively more improbable. In fact, for e sufficiently large, the false vacuum may 
be below percolation threshold and only the true vacuum will percolate. This situation 
was analysed previously in the context of the formation of NTSs 8 but will clearly not be 
of interest here; we must insure that both vacua percolate even though they will have 
different populations. This will effectively put an absolute upper bound on e which is 
obtained below. 

As long as the system is in equilibrium, the relative population of the two vacua is given 
by the Boltzman formula, p+ /p_ = exp(— AF/T), where A F = A x is the difference in 
free energy between the two minima. As explained above, these populations will be frozen 
as the temperature drops below Tq. Since Tq ~V^x Um, for T < Tq we obtain, 

P+/P- = exp(— A/Ifjtf) sf exp(— 4e/A), (14) 

where the last result is obtained by taking A cs egq and Um ^ (A/4)<7q. Notice that we 
are using the expression for Um in the limit where e — > 0. This approximation is valid as 
long as Um A, which is reasonable to expect for small asymmetry. We can now obtain 
the upper limit on e that results in percolation of both vacua. If one considers a simple 
cubic lattice (which we will, as a working model), p c = 0.31. 10 Thus, solving Eq. (14) for 
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p+ > 0.31 one obtains 

e < 0.2A. (15) 

Other values for p c would just change slightly the numerical coefficient in front of A without 
changing the linear behaviour of e = e(A). With this upper value satisfied, even though 
the vol um e occupied by the true vacuum will be larger than that occupied by the false 
vacuum, at the completion of the phase transition space will be permeated by an infinite 
domain wall separating regions of plus and minus vacua. 

We conclude this section by calculating the energy density in walls at formation, i.e., 
at T = Tq- In the last section we obtained the average energy density in walls to be 
pjy ~ p/(L). We here give a more precise derivation. The wall energy density, p\y, is 
proportional to the surface tension p. The proportionality factor is just the ratio between 
the area in walls over the volume, A wa u/V. If we divide a space of volume V into N 3 cubic 
cells of volume £ 3 each, the probability of having a wall between two cells is given by twice 
the probability of having a plus-cell p+ followed by a minus-cell 1 — p+ = p~, i.e. 2p(l — p). 
Clearly, this is also the ratio between the area of the walls in the volume V to the area of all 
cell-boundaries in V, ^wallMtotal* Ih order to obtain A wa u/V we must calculate the ratio 
A otal/^j i- e - the ratio between the total area of cell boundaries over the volume. As we 
have N 3 cells of volume £ 3 each, there will be N — 1 internal boundaries between any two 
cells perpendicular to the x-axis, each with area JV 2 £ 2 , the same happening for the y and z- 
axis. Thus, the total area in cell boundaries (i.e. boundaries between two cells, where walls 
may appear if the vacua on both sides are different) is Aotal = 3(AT — 1 )AT 2 £ 2 — 3N 3 £ 2 , 
where in the last expression the large N limit was assumed. The desired ratio is then 
AmJV = 3N 3 ( 2 /(N(f = 3/e. Thus, 

-^wall ^wall -^total _ ^ _\ fm\ 

(16) 

Using the value of p obtained in Eq. (6), the approximate value for the correlation length 
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at Tq, £{Tg) ~ (Actq) - 1 , and Eq. (16), we find for the density in walls at formation 

Pw{Tg) - 4\/^p(l - p)A 3 / 2 <Tq. 

Thus, at r<3, the ratio of energy density in walls and radiation is [recall that pp = g*Tg 
<7*<7q( A + l/4)“ 2 , where <7* is the number of relativistic degrees of freedom at Tq\, 

Or r " + V<) 3/2 - (18) 

' PR ' T=Ta g* g * 

As is well known (see Refs. 4-7, 11), the walls do not dominate the energy density of the 
Universe at formation. 

4. Dynamics of Walls I: Convoluted Structures Inside the Horizon 

As the walls are formed at Tq , space will be filled by the complicated wall structure 
described in the last section. Typically, the average wall separation ( L ), and the average 
curvature radius of the wall, R, will be of the same order as the correlation length £(T <y) c* 
(Aao) -1 . The motion of the walls will be determined by two forces. 

The first force is from the surface tension. Due to curvature the walls will feel a 
surface pressure that acts to straighten them. So long as the wall curvature is smaller 
than the horizon size, i.e., R < A jj cs Mpig^^ctQ 2 where Mpi is the Planck mass, 
this straightening will be rapid compared to the age of the Universe at that time. The 
surface pressure is given by pp ~ fj,/R(t). At formation, pp ~ A 3 / 2 Oq. If this was the 
only force acting on walls, 4 it would induce oscillations of frequency v ~ £ , quickly 

accelerating the walls to relativistic speeds. Of course, the motion of walls will be damped 
by radiation of particles and by the interactions of walls with each other and with the 
medium. 5 As they evolve, the walls will become progressively flatter and slower up to 
the point where their curvature is of the same order as the horizon scale. Note that this 
damping mechanism is occuring on the scale of microphysics (we assume that oq <C Mpi) 
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and that the flattening of the walls will take a negligible time compared to the cosmological 
time scale. (For the readers convenience we include in Table 1 a complete list of all 
time scales used in this paper.) It is then easy to estimate the time when the walls can 
dominate the energy density of the Universe causing a power-law inflation (Refs. 4-6, 
11), p\y = A x l 2 o%/t > PR ~ for t > tjnf a* A -1 / 2 ^ 3 . Let us call the 

time at which power-law expansion starts, and tg the time of wall formation, for which 
T = Tg. Thus their ratio is simply t^/tg =* {gV 2 / A 1 / 2 )(Mp i /ao). This result was the 
root of the domain wall problem; the walls would eventually dominate the energy density 
of the Universe, causing a power-law inflation and conflict with astrophysical observations 
as explained in the introduction. Walls had to disappear before t' m f, or possibly some 
time not much larger than - 11 In the present scenario though, the walls will also suffer 
the influence of the volume pressure due to the asymmetry. The dynamics will be quite 
different. 

The second force is a volume pressure on the walls due to the energy difference between 
the two vacua, py ~ eo o- This pressure will accelerate the walls against the false vacuum 
regions, rapidly converting false vacuum into true vacuum. The energy released in the 
process fuels the wall motion. The process is very similar to the propagation of relativistic 
detonation waves analysed by Steinhardt. 9 The wall will reach some terminal velocity io e q 
that will be dependent on the type of matter on both sides of the wall. (This is the 
wall velocity discussed in Section 2. The subscript eq was included now for clarity.) For 
example, the wall may be converting regions where a gauge symmetry is unbroken into a 
broken symmetry phase by giving mass to a certain particle, or the wall may be moving 
between relativistic fluids with different chemistry. The dynamics of the walls will be very 
sensitive to both the pressure and the tension forces; the exact treatment being out of the 
scope of this work. 

The dynamics of the wall system will be examined in the various possible regimes that 
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will depend upon when and which of the two forces dominate. 

It is clear that the initial dynamics of the walls will depend on the balance between the 
two forces at formation. From the above, the ratio between the two pressures at t = tQ is 

(?£.) ~ a- 3 /2. (is) 

\prh-ta 

There will be two possible regimes. For py > pt e > A 3 / 2 ), the volume pressure will 
dominate the motion of the walls from the moment they are formed; for py < PT 0- e -> 
e < A 3 / 2 ), the tension will dominate the motion of the walls initially. Let us call these two 
regimes A and B respectively. Case B can be further divided in two subcases, depending on 
whether the volume pressure acts on the walls before or after the tension has straightened 
them to the horizon scale. We call these two cases B.l and B.2 respectively. Case B.2 is 
quite involved and will be treated separately in the next section. We proceed here to treat 
cases A and B.l, when the wall is still convoluted inside the horizon when py dominates 
PT- In Table 2 we present a brief summary of the four cosmological scenarios examined in 
this paper. 

In addition to these conditions determining the initial behaviour of the walls, we must 
also ensure that both domains percolate by satisfying the inequality obtained in Eq. (15). 
(Otherwise the Universe would consist of a true vacuum sea with isolated bags of false 
vacuum that would, in the absence of stabilizing pressures, quickly disappear due to the 
tension force, as explained in Ref. 8.) In Fig. 2 (that should be interpreted together 
with Table 2), we show schematically the range for each of the two regimes with the 
percolation requirement satisfied. Notice that the two curves cross when A 3 / 2 = 0.2A. 
The fact that case A is given by a lower bound for e imposes a strong constraint on the 
range of allowed parameter space; A < 0.04 in this case. The walls will be immediately 
accelerated against the false vacuum regions by the volume pressure. The acceleration is 
roughly, ay ~ py /p = ecro/A 1 / 2 . As the walls are separated initially by ( L ) ~ £{Tq), the 
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time necessary for them to pinch off is ~ (A 1 / 2 er) 1 / 2 cr 0 - They will disappear before 

any noticeable effect can occur on the cosmological time scale. For example, if A = 10 -2 , 
tpinch/^G ~ 10 3 ao !Mp\. There is no need to invoke Steinhardt’s analysis for detonation 
waves because the time required for the shock front to acquire equilibrium, £ e q ~ Weq/ay, 
making his analysis applicable, is larger than tp^ch by a factor of order e -1 / 2 . 

We can now move on to case B. The conditions on the asymmetry are e < A 3 / 2 and 
e < 0.2A. Thus, contrary to case A, only upper bounds are imposed on the asymmetry; as 
we will show, interesting scenarios are possible for small values of e. In case B.l, the volume 
pressure becomes important before tension can straighten walls in a horizon scale but not 
immediately after formation, as in case A. As py becomes dynamically effective, there will 
still be an entangled structure within a horizon size, although the average curvature radius 
and wall separation will now be bigger than at t G by an amount that scales with e _1 . 
The condition that the volume pressure becomes effective before the walls are straight on 
a horizon scale imposes a lower bound for e: e > A 1 / 2 cr^" 1 ^(f) -1 . For R(t) ~ A#, we get 

£ > (As - )I/2 te) • < 2 °) 

In Fig. 2 we show graphically the range of parameters that satisfy the three inequalities for 
case B.l. The motion of the walls is initially dominated by the tension force that will make 
them straight on scales up to R{t e ) = A^/eco < Ajy, when volume pressure dominates 
and starts accelerating walls toward the false vacuum regions. The disappearance of the 
walls in this case is very similar to case A, the difference being only on the scales. At 
t e the walls will be separated on average by R(t e ). The walls will disappear after a time 
*piLh — Using (20) we obtain that t pinch < M Pl g* l ! 2 o q 2 = 0(t G ), where t G is 

the age of the Universe at formation. So the pinching occurs before the walls can have any 
major role in the cosmological evolution. Steinhart’s analysis is still not relevant, although 
now t e q is bigger than in case A; f e q ~ A 1 / 2 /eo'o ^ * p inch- The walls start pinching as 


18 



IF? I 1 ,a l ® ^ S'* .■■■•> 

OF tOljit '/'i* k .-*v 

they reach terminal velocity, with the time scale for both processes being smaller than the 
cosmological time scale. 

5. Dynamics of Walls II: Structures Curved on a Horizon Scale 

In the last section we looked at the dynamics of walls in the case where the volume pres- 
sure coming from the asymmetry dominates the wall motion before the tension force can 
straighten the walls on a horizon scale. We have seen that the asymmetry causes the dis- 
appearance of the walls before they can play any major role in the cosmological evolution. 
In the present section we will consider the case when the asymmetry is sufficiently small 
in order to have walls with average separation and curvature comparable to the horizon 
size before the volume pressure starts to become effective. 

This situation can still be subdivided into two cases (which are called cases B.2.1 and 
B.2.2 in Table 2 and in Fig. 2). One possibility, case B.2.1, is to have the volume pressure 
accelerating the walls before they can dominate the energy density of the Universe. This 
situation is the one explored earlier in the literature (see, e.g., Refs. 4-7). The other 
possibility, case B.2.2, is to have the walls dominating the energy density of the Universe 
before the volume pressure can act. Although the possibility of a wall-dominated Universe 
is perfectly reasonable, it is not clear to us what the evolution of such Universe will be 
like. We will present a rather naive approach to this case, being careful to point out the 
serious difficulties involved. 

In order to have walls with curvature of the same order as the horizon we need simply 
to invert the inequality in Eq. (20), 

E< < A »-> 1 / 2 te)’ ' 21 > 

with the other two conditions on e the same as in case B.l (again, shown in Fig. 2). Once 
these three conditions are satisfied, the volume pressure will only become dynamically 
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effective after the tension has straightened walls on scales up to \jj. We still need to 
distinguish between the two possibilities discussed in the previous paragraph by imposing 
suitable bounds on e. For this we simply compare t e , the time when the vacuum pressure 
dominates wall motion, to the time when the walls dominate the energy density of 
the Universe. From the results in the previous section we get 

*inf e f M Pl \ 2 

t e A \ oq ) 

Thus, in order to avoid a wall-dominated Universe we must obey a lower bound for the 
asymmetry, e > \(oo/Mpi) 2 . This lower bound has been found before by Vilenkin. 6 It 
must, however, be supplemented by the other two upper bounds on e given in Eqs. (15) 
and (19) for this scenario to work. Assuming that we live in the true vacuum, the walls will 
be moving away towards the false vacuum with relativistic speeds that can be calculated 
using Steinhart’s analysis. 9 For us, it is important to note that the terminal velocity will 
be in most cases l/-\/3 < tu e q < 1, the exact value depending on the “chemistry” of the two 
phases. Of course, as the Universe expands and the temperature drops, wall velocities will 
asymptotically approach the speed of light. Such walls should not cause any cosmological 
damage. 

The second possibility is to have the walls dominate the energy density of the Universe 
before the volume pressure turns on. From Eq. (22) we can see that this case obtains for 
exceedingly small asymmetry (unless, of course, <tq ~ Mp/), that is, for almost degenerate 
vacua. As has been shown in Ref. 4 such Universe would have a power-law inflation with 
the scale factor evolving like iZ(t) ~ t 2 . We have shown in Section 2 that this is indeed true 
in the limit of having a very non-relativistic gas of walls in a box. We question, however, 
if our analysis can be applied without contradictions to the situation with one wall per 
horizon volume. It seems that the more adequate approach to this situation would be 
to solve Einstein equations for one wall in a Robertson- Walker Universe. To the best of 
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our knowledge, this has not been attempted so far. The possibility of a wall-dominated 
inflation remains an open question. Nevertheless, let us, for the sake of an argument, 
naively assume that wall-dominated inflation is possible once there is one wall per horizon 
and see what we can learn about the role of an approximate discrete symmetry in such 
Universe. 


Assuming that the upper bound for e < X(ao/Mpi) 2 is satisfied, the walls will be 
able to dominate the energy density of the Universe before the volume pressure turns 
on. Assuming further that the t 2 behavior for the scale factor is correct, we can ask how 
small e has to be in order to have sufficient inflation to solve the horizon and flatness 
problems of standard cosmology. 12 In this naive picture, inflation would end once the 
walls are accelerated by the volume pressure to velocities larger than tu eq = l/\/3> since 
for u>eq > l/\/3 the scale factor will evolve slower than the horizon, as shown in Section 2. 
As the walls are accelerated with ay ~ eao/A 1 / 2 it is easy to show that £ e ndinf> the time 
when inflation ends, and t € are of the same order in comparison to the cosmological time 
scale. A rough estimate for the maximum (the exact calculation is quite involved since the 
exponent for the time evolution of the scale factor is itself a function of time) amount of 
inflation is then, 


-^foendinf) < ( ^endinf \ A 2 / Op \ 

*(*inf) “ \ *inf / ~ * 2 \MPI) 


(23) 


In order to solve the horizon problem and make the walls disappear, we need 
•ft(*en dinf )/-fl(*inf) ^ 10 26 . 12 Thus, the upper bound on e is 


e < 10 -13 A 



(24) 


A successful wall-dominated inflation would require, in this grossly simplified calculation, 
an extremely small asymmetry! Of course, even assuming that this analysis is qualitatively 
correct, only through simulations of wall motion including damping from the medium and 
the Hubble flow could a more explicit picture emerge. 
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6. Concluding Remarks 

The role of topological defects in cosmology was first studied in the work of Zel’dovich, 
Kobzarev and Okun. 4 This work, and the extensive literature that it has inspired, has 
influenced the direction of model building in particle physics, in particular when discrete 
symmetries were involved, as for example in axion models and CP violation models. The 
dr a.ma. tic results of Ref. 4, that particle theories with exact discrete symmetries are ruled 
out by cosmology, can be relaxed, as suggested in that very work, in the case the discrete 
symmetry is only approximate. As we hoped to have shown in this work, this is indeed 
the case. 

We studied in some detail the cosmological consequences of having an approximate 
discrete symmetry spontaneously broken in the early Universe. We have shown that many 
possible scenarios are possible, depending on the value of the asymmetry parameter. Ba- 
sically, after domain walls are formed in a primordial phase transition, their dynamics will 
be determined by the interplay between the tension force due to the curvature of the walls 
and the volume force coming from the volume energy stored in the false vacuum due to 
the asymmetry. 

We divided the possible cases according to which of the two forces dominates first. 
Thus case A has the volume force dominant from the walls appearance while case B has 
the tension force dominating the dynamics initially. In case A the walls disappear in a 
time much smaller than the cosmological time scale. Case B was divided into two subcases, 
depending on whether the tension force has acted to straighten the walls up to the scale 
of the horizon at the time when the volume pressure starts acting. In case B.l the volume 
force acts before this can happen, and the walls disappear as they begin to be of any 
cosmological relevance. In case B.2, there will be essentially one wall per horizon before 
the volume pressure can act. This case can be further divided into two subcases, according 
to whether the energy density of the Universe is wall dominated or not. In the latter 
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case, case B.2.1, the volume pressure will accelerate the wall to relativistic speeds and 
the Universe will be radiation dominated. In the first case, case B.2.2, if one accepts a 
naive approach to the problem of a wall-dominated Universe, we found that there can 
be an inflationary epoch that ends as the volume pressure starts accelerating the walls 
towards the false vacuum regions. Each of these scenarios obtains for certain bounds on 
the asymmetry that we evaluated in terms of the parameters involved in the simple model 
considered here. 
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Table Captions 

Table 1. A time table for the different time scales defined in the text. 

Table 2. A table of scenarios considered in the text. Recall that for formation of 
infinite walls, e < 0.2A 
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Figure Captions 


Fig 1. The zero temperature potential for non-degenerate vacua, e, the asymmetry 
parameter, is in general function of temperature and of different couplings in the model. 


Fig 2. Schematic display of all four possible scenarios for different values of the asym- 
metry parameter e. In case A the walls disappear in a time scale much smaller than the 
cosmological time scale, tQ. In case B.l the walls disappear in a time scale comparable to 
tQ. In case B.2.1 there can be one straight wall in a horizon volume but the walls disappear 
before dominating the energy density of the Universe. In case B.2.2 the walls dominate 
the energy density, causing a power law inflation until they disappear. 
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Time scale 

Definition 

Time when... 

*G 

“ M Pl 5 r 1 / 2 c 0 - 2 

walls are formed. 

^inf 


walls dominate energy-density 
of Universe. 

Spinet 

o' (A 1 / 2 e) -1 / 2 a^' 1 

walls pinch off. 

teq 

a ^cqA 1 / 2 (£rao) -1 

walls reach terminal velocity w e q . 

t e 

a A 1 / 2 (ecr 0 )— 1 

vacuum pressure dominates wall motion. 

*endinf 

a A 1 /2(ecr 0 )-l 

wall dominated inflation ends 
(valid only for e < A(oo/Mpi) 2 ). 
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TABLE 2 


Case 

Scenario 

Bounds On Asymmetry 

A 

Volume pressure dominates 

initially. Walls 
disappear very fast 

( U ct G ). 

£ > A»/2 

B.l 

Tension force dominates 

initially but cannot straighten 
walls up to horizon scale 

(te < t G ). 

A 3 / 2 > e > (Aj7*) 1 / 2 (a 0 /Afpi) 

B.2.1 

Walls are curved on 
a horizon scale but cannot 
dominate energy-density 
(*e < ^inf)* 

e < A 3 / 2 ; e: < (Ap*) 1 / 2 (o 0 /M P i) ; 
e > X(oo/Mpi) 2 

B.2.2 

Walls dominate energy-density 

causing power law 

inflation before disappearing 

e < A 3 / 2 ; e < {Xg*) l l 2 {o Q lMp\) ; 
e < A(<7o/Afpi) 2 
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